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Abstract 

Light pseudoscalars known as axion like particles (ALPs) may be behind physical phenomena like the 
Universe transparency to ultra-energetic photons, the soft y-ray excess from the Coma cluster, and the 
3.5 keV line. We explore the connection of these particles with the inverse seesaw (ISS) mechanism for 
neutrino mass generation. We propose a very restrictive setting where the scalar field hosting the ALP is 
also responsible for generating the ISS mass scales through its vacuum expectation value on gravity induced 
nonrenormalizable operators. A discrete gauge symmetry protects the theory from the appearance of overly 
strong gravitational effects and discrete anomaly cancellation imposes strong constraints on the order of the 
group. The anomalous U(l) symmetry leading to the ALP is an extended lepton number and the protective 
discrete symmetry can be always chosen as a subgroup of a combination of the lepton number and the 
baryon number. 
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I. INTRODUCTION 


Pseudo Nambu-Goldstone bosons of very low masses, interacting effectively with the electro¬ 
magnetic field, are generally predicted in Standard Model (SM) extensions containing approximate 
continuous global symmetries which are spontaneously broken. A distinctive example of this type 
of particles is the axion, which arises when the strong CP problem is solved through the Peccei- 
Quinn mechanism m- Generically, any light pseudo Nambu-Goldstone boson whose coupling to 
photons is similar to that of the axion has been denoted as an axion like particle (ALP). Several 
experiments are in search of the peculiar effect of photon o ALP oscillations and these searches 
guide the development of theories containing these particles [1H6]. 

The symmetry associated to the ALP is usually taken to be spontaneously broken at a very 
high energy scale, and we assume here that this breakdown occurs through the vacuum expectation 
value (vev) of a scalar field singlet under the SM symmetry group. The ALP decay constant, which 
controls the feeble ALP couplings to other SM particles, can be identified to this vev, up to order 
one coefficients. 

In this work we investigate settings where the same scalar singlet hosting the ALP also gives rise 
to the mass scales involved in the inverse seesaw (ISS) mechanism for neutrino mass generation, 
a well-known mechanism that does not require too heavy neutral fermions [7]. There are two 
scales beyond the SM in the ISS mechanism. One of them is directly related to the lepton number 
breakdown. In our proposal the usual lepton number is embedded in an anomalous U(l)x symmetry 
associated to the ALP. Such a symmetry is an accidental one in the sense that it automatically 
arises from the imposition of another symmetry considered to be more fundamental. Breaking of the 
accidental U(l)x symmetry is expected from gravitational interactions through nonrenormalizable 
operators, which also generate the ALP mass. 

Therefore, we also deal with the problem of stabilizing both the mass scales in the ISS mechanism 
and the ALP mass in face of gravitational interactions. Nonrenormalizable operators that would 
bring too large mass corrections can be avoided by assuming discrete Z^v symmetries, which are 
remnants of gauge symmetries valid at very high energies [8]. The choices for the Z^r symmetries 
are greatly reduced by the conditions they need to satisfy in order to be free from anomalies mm- 
For example, there are only a handful of interesting anomaly free discrete gauge Ztv symmetries 
for the MSSM P ITT]. 

Discrete symmetries have already been used for suppressing dangerous operators that prevents 
the solution of the strong GP problem in different models containing the axion |12H21j . Such 
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symmetries are shown to be crucial in multi-ALP models, where very low mass ALPs need effective 
protection against gravitational interactions [ 20 ) . 

At the same time that the U(l)x symmetry breaking scale leads to the correct mass scales for the 
ISS mechanism, we look for models where the ALP coupling with photons and the ALP mass have 
values that allow the explanation of three hinted astrophysical phenomena: the anomalous Universe 
transparency for very energetic 7 -ray [MSB, the soft X-ray excess from the Coma cluster [32H36j , 
and the X-ray line at 3.5 keV |37H40j . These phenomena have already motivated the development 
of general multi-ALP models, containing an axion dark matter candidate whose decay constant is 
associated to the high energy scale entering in the canonical seesaw mechanism [ 2 U| . 

The first astrophysical hint, the anomalous Universe transparency, follows from observations 
of active galactic nuclei (AGN) |24H27[ 131| . An ALP coupling with the photon and with mass 
within some specific interval could provide an explanation to this phenomena through 7 -ray o 
ALP oscillations in the magnetic fields of AGNs, and in the intergalactic medium |23l I28ll30j (see, 
however. Ref. |1B ^ recent analysis of the cosmic transparency hint). 

Second, the observed excess of soft X-ray coming from the Coma cluster |32j could be explained 
assuming a cosmic ALP background radiation, originated from the decay of heavy moduli fields 
and corresponding to a fractional number of extra neutrinos |33l [34] . The cluster magnetic field 
would make possible the conversion of the ALP into the observed X-rays |35[ (36] . 

At last, the 3.5 keV line has been reported from observations, using the XMM-Newton satellite 
data, of the stacked spectrum of galaxy clusters and in the Perseus cluster |3Z|, and also in the 
Andromeda galaxy |3B- A possible interpretation for this could be the two photon decay of a 
dark matter ALP with mass equal to 7.1 keV |39l 140] (another possibility could be a CP even 
scalar [42] . or a specific majoron[43|). Even if a 7.1 keV dark matter does not correspond to an 
ALP, it could decay into ALPs that in turn decay into photons [44] . It has to be said that a 
study with the Chandra data of X-ray of the Milky Way did not show a conclusive evidence for 
the 3.5 keV line [45|, and that other interpretations for the 3.5 keV line in terms of some specific 
Potassium and Chlorine lines were also suggested |46] . On the other hand it is argued in m that 
the interpretation of dark matter decay as the origin of the 3.5 keV line is consistent with the 
XXM-Newton dataset of the Milk Way center. Although there is some debate on the origin of the 
3.5 keV line signal we shall assume it is due to an ALP decay. 

We also motivate our work with the new generation of proposed experiments which are projected 
to probe some regions of the parameter space of the ALP coupling with photons and its mass. 
Among these experiments, we can mention: the ALPS-II [48], the helioscope lAXO [49], and the 
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observatories PIXIE [50] and PRISM [5T]. Most of the models we propose here are within the 
prospected search range of these experiments. In Figure]^ it is shown the regions in the parameter 
space to be tested by these experiments, as well as the ones allowing for explanation of the hinted 
astrophysical phenomena. 

The outline of the paper is the following: in Section [ll] we present the general setting that 
relates the physics of an ALP, its astrophysical motivations and the generation of the ISS scales. In 
subsequent Section [II A we analyze the general symmetry properties of the models and establish 


necessary conditions for interesting models. Then in sections IIB and IIC we show, respectively. 


models with one and two ALPs. Finally, We conclude in Section III 


II. ALP AND THE INVERSE SEESAW MECHANISM 


We start by showing the main elements that need to be considered in our constructions con¬ 
taining just one complex scalar held whose vev generates the energy scales involved in the ISS 
mechanism, and which are assumed to be associated with new physics beyond the SM. Such scales 
are taken as being proportional to a scalar held vev times a suppression factor, composed by this 
vev divided by the Planck scale and raised to some power. The complex scalar held hosts an ALP 
which, through its effective interaction with the electromagnetic held, is going to provide expla¬ 
nation for astrophysical phenomena like the soft X-ray excess and the Universe transparency. For 
this, the ALP needs to have its mass protected from dangerous effective operators due to gravita¬ 
tional interactions. In order to obtain the natural ISS mechanism scales and the appropriate mass 
for the ALP, we look for suitable discrete symmetries over the helds. 

In the ISS mechanism [7], two extra sets of neutral fermionic singlet helds, Niji and SiR, i = 
1,2,3, are taken into account in addition to the SM neutrino helds UiR. It is assumed that after 
spontaneous symmetry breaking, a mass Lagrangian is generated containing the following terms 

D Nr m£) i/R + Sr M Nfi + 2^^ ^ 


with the 3x3 Dirac mass matrices nriR, M, and the Majorana mass matrix which without loss 
of generality can be taken diagonal. The mass matrix texture arising from Eq. Q, with the basis 
choice [h'L , N^ , is 


M = 


0 mj) 0 

niD 0 M"*" 

0 M jj 


( 2 ) 
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It was observed by Mohapatra and Valle [7] that a mass matrix of the form in Eq. Q may lead 
to three active neutrinos with masses at the sub-eV scale, without invoking very large entries in 
the matrix M. For example, masses at the sub-eV scale for the active neutrinos can be obtained 
with the entries of m/j, M, and /r of order 10 GeV, 1 TeV, and 1 keV, respectively. Specifically, the 
lepton number is only broken by a small scale set by /r, which is the inverse of what is assumed 
in the canonical seesaw mechanism where lepton number is broken by a very large right-handed 
neutrino scale. Taking a matrix expansion in powers of M~^, block diagonalization of M leads to 
the approximate mass matrix for the three light active neutrinos 

-1 


^ mjj M tud , 


( 3 ) 


and a 6 X 6 matrix 


Mr 


0 

M jji 


( 4 ) 


related to six neutrinos. These last ones are supposedly heavier than the active neutrinos, with 
masses at the scale of M, and are quasi-degenerate (pseudo-Dirac nature) if the entries of // are 
small compared to the ones in M. If the number of Sir fields were greater than the number of 
NiR fields, one or more neutrino states with masses at the ^ scale would arise, and they could 
also contribute as dark matter [52] (another possibility for keV DM within the ISS mechanism is 
given in Ref. [53]). The mixing between the heavy neutrinos and active neutrinos is approximately 
given by e = m£) M~^ and unitarity violation effects are typically of the order e^. General aspects 
of the ISS mechanism concerning the neutrino mixing and violation of unitarity were developed 
in Refs. [5ll[55|- Generically, at the percent level is not excluded experimentally, but may be 
within the reach of future experiments probing lepton flavor violating transitions m and direct 
production of heavy states at colliders [56] . 

The scales involved in M and /r are supposed to arise from new physics beyond the SM. In 
particular, the fi term in Eq. 0 breaks the lepton number symmetry explicitly. From the point 
of view of naturalness it is reasonable that the nonvanishing entries of /r be associated with a 
small effective energy scale compared to the electroweak scale, Uw = 246 GeV. In the limit /r —)• 0 
lepton number conservation is restored increasing the set of symmetries. Thus, the entries of /U 
are expected to be small in comparison with the mass scales of the SM, which contain the lepton 
number as a global automatic symmetry. 

In our approach, the parameters /r and M are gravity induced and result from the very high 
vev of the complex scalar field times suppression factors containing the Planck scale, with the 
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parameters in mu proportional to v^. This contrast with proposals where the typical energy scale 
in M is due to a new theory with spontaneous symmetry breaking at the TeV scale [57H59| . 

The complex scalar field is a singlet under the SM gauge group and has a vev denoted as 
(cj) = VajV^, with 

109GeV<u, < IQi^GeV, (5) 


defining the intermediate scale range. This leads to a photon-ALP coupling with value required 
to explain astrophysical phenomena, with the ALP detectable by future experiments |20l l60] . We 
parameterize the scalar singlet as 


. a{x) 

= — [va + p{x)]e^ 


( 6 ) 


were a{x) is the ALP field. The radial field p{x) gets a mass at the scale and we assume it 
decouples from the low energy effective theory. In the models presented here a carries charge 
of a global U(l)x chiral symmetry which is explicitly broken by the gravitational interactions in 
such a way that, after spontaneous symmetry breaking, the ALP gets a small mass. The U(l)x 
symmetry is taken as accidental meaning that it results from one or more imposed gauge discrete 
ZAT symmetries - not broken by gravitational interactions - restricting the main interactions of the 
neutral fermion fields with the scalar fields being 


^ Nr y H^L + ^p-i V Nr + ^ Sr C Sr + H.c. , 


(7) 


where y, rj are complex 3x3 matrices, and C is a symmetric 3x3 matrix. Li and H are the 
leptons and Higgs SU(2 )a doublet fields, respectively, with H = iT 2 H*. The complex conjugate 
field a* can be equally considered in the third term, instead of a, while we conventionally define the 
scalar present in the second term to be a. We use the reduced Planck scale Mpi = 2.4 x 10^® GeV 
for the gravitational scale. The vev of the Higgs doublet field is (H) = (0, Uw/\/2)^. We will 
see that U(1 )a: is directly related to an extended lepton number and thus, in our approach, the 
smallness of p follows from its explicit but small breaking due to gravity (1/Mpi suppression) and 
its spontaneous breaking at the scale Va- 

With the vev of the scalar fields in Eq. the effective Lagrangian in Eq. 0 is obtained with 
the mass matrices 


ruD = yNN , M = r]- 


'V2' ''2fMSr'’ 


p = c 


^P1 


2iM|“^ 


( 8 ) 


niD is naturally at the 100 GeV scale without requiring the entries of y to be fine tuned. The non¬ 
vanishing entries of rj and are all of order one, under the assumption that the nonrenormalizable 
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interactions in Eq. 0 are exclusively due to gravitational interactions, whose universal coupling 
is 1/Mpi. Thus, we look for values of in which the mass scale function F{k) = v^/22 is 

assumed to have values F{p) = 0.1 - lOTeV and F(q) = 0.1 - lOkeV, for p and q integers. In 
figure curves for p and q are shown as functions of v^. Within the range in Eq. 0 , we can see 
that only 

p = 2,3 and q = 3,4,5, (9) 


can generate the appropriate scales for the ISS mechanism. Moreover, if only one vev accounts for 
both M and p, it can be seen that there are only two sets of solutions: 


(p, q) = (2, 3) for Va ~ 2.2 x 10^° GeV — 5.5 x 10^° GeV, 
(P) q) = (3, 5) for Va ^ 2.8 x 10^^ GeV — 5.5 x 10^^ GeV. 
However, we should keep in mind that the scale for p is more flexible than M. 



( 10 ) 


FIG. 1: p (purple) and q (red) as a function of for F{p) = 0.1 - lOTeV and F{q) = 0.1 - lOkeV, with 
F{k) = uJ/ 25The regions of intersection between yellow and blue bands indicate common values for 
Va, possible for the set of integers {p, q) = (2, 3) or {p, q) = (3, 5). 


We observe that for both sets of values of (p, q) in Eq. (10) the active neutrinos masses are 


independent of the Planck mass, at leading order, and it happens whenever q = 2p — 1. This feature 
would not be possible if, e.g., neutrinos masses were generated by type-I seesaw with heavy masses 
induced by gravity. In fact, Eq. 0 leads, in face of Eq. ([^, to 




y^v ^C{v 


-1 


V2i 


( 11 ) 
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This formula shows that the active neutrinos masses have a suppression factor v^jvo- in relation 
to the electroweak scale, v^. Such an explanation for having small neutrinos masses resembles 
the canonical seesaw mechanism and have been observed in other models implementing the ISS 
mechanism m- 

Now we show that besides having its value constrained to be within a range that allows active 
neutrinos to have masses at the sub-eV scale, can also furnish an ALP-photon coupling as 
required for explaining the mentioned astrophysical phenomena. It is supposed here that such 
ALP-photon coupling arises effectively by means of the a field interaction with a new colorless 
fermion, E, which carries one unit of electric charge and is a singlet under SU(2)/^ group. Along 
with (T, NiR, and Sir, we assnme that E also carries charge of U(l)js:- Under this symmetry these 
fields transform as 


NiR —)> E^^^NiR , SiR —)> e^^^^SRi, 
El ^ E^^lPEl , Er ^ E^^rPEr, 


( 12 ) 


with charges A^, normalized such that X^j = 1. 

The effective Lagrangian describing the interactions of the ALP with photons is 


— 


\d,,ad^a — a? 


9a'y 

X 




(13) 


where is the electromagnetic field strength and its dual. The ALP-photon 

coupling constant, ga'y, is given by 


9a'y 


^ ^a'y 
2 tT Vfj ’ 


(14) 


where a ^ 1/137, and the anomaly coefficient Ca^y depends on the U(l)x and electric charges of 

(i) 

the fermionic fields 'ipi, ^ and Cem, respectively, according to 


Ca, = 2Y,{X^,-X^^){Cii^f. 


(15) 




Throughout this work only the field E is chiral under U(l)x “ the left- and right-handed components 


of E have different U(l)x charges - so that Eq. (15) reduces to Ca^y = 2{Xe^ — Ag^). The effective 


coupling in Eq. (14) can be obtained through a rotation of the fermionic fields ip —)■ ip which 


a(a;) 


does not leave the integration measure invariant, meaning that the U(l)jv symmetry is anomalous 
(for details see Ref. [HI])- With such a rotation it turns out that the ALP has only derivative 





couplings with the fermions. We omit the interactions of the ALP with fermions since these effects 
are outside the scope of this work. 


A nonzero value for nia in Eq. (13) must be generated by an explicit breaking of U(l)x, 
characterizing the ALP as a pseudo Nambu-Goldstone boson. We follow the assnmption that gra¬ 
vitational interactions do not respect any global continuous symmetry and that U(l)x is explicitly 
broken by nonrenormalizable operators suppressed by Mpi. However, as argued by Krauss and 
Wilczek [ 8 ], discrete symmetries like Zjv, which are leftover of gauge symmetries, are expected 
to be respected by gravitational interactions ^ and, therefore, they can prevent the presence of 
unwanted nonrenormalizable operators of lower dimensions. Thus, a Zjv preserving operator of 
some high dimension D necessarily breaks U(l)x, 


.if D 


M, 


^ rCr'^-l-H.c. 


D-4^ 


p>4) 


(16) 


PI 


where g = |g'|e*'^, with \g\ of order one; we assume the operator in Eq. (16) is the one with lowest 
dimension with such a property. In that case, at leading order, a potential for the ALP is generated 


V{a) 


|g| Vg 
2 (\/2Mpi)^-4 


cos 


D— + 6 

Vrj 


(17) 


This furnishes a mass to the ALP which can be very light for a sufficiently high D, 

(18) 

Intervals for the ALP parameters {ga'y, ma) which can explain the anomalous Universe trans¬ 
parency for very energetic 7 -ray [22H3U| . the soft X-ray excess from the Coma cluster 133 ES]: 
and the X-ray line at 3.5 keV [3711381130], are shown in Figure It can be seen that there is a 
region where a set of parameters could explain both the anomalous Universe transparency and the 
soft X-ray excess from the Coma cluster. That region corresponds to ga'y ~ 10“^^ — 10“^^ CeV“^, 
which implies an ALP scale at the range 


m„ 






PI 


^ pa 10® — 10^ CeV . 


(19) 


07 


with tUq < 10 ^^eV. This requires that the U(l)x breaking operators in Eq. (16) should have 


dimensions of at least D = 11 for Va = 10^ CeV, and D = 12 for Va = lO^'^ CeV 


^ The argument in Ref. [S] is essentially that gravitational interactions must respect local symmetries and also any 
residual Zjv symmetry left in the effective theory after spontaneous breaking. 
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On the other hand, the region of parameters allowed for explaining the X-ray line at 3.5 keV is 
disconnected from the previous region. Thus, if all these hinted phenomena are due to ALPs at least 
two different species of them are needed to exist. In order to explain the 3.5 keV X-ray line through 


a decay of an ALP with mass 7.1 keV, the dimensionality of the operator in Eq. (16) inducing such 


a mass depends on the scale Va. The values for the ALP scale and the operator dimension D 
inducing the correct mass with coupling in the range Qa-yfCa-y ~ 10“^^ — 10“^^ GeV“^ are shown 
in Table For simplicity we take <7 = 1. We see that the U(l)x breaking operator needs to be of 
dimension 7 or larger. 

Large portions of the ALP parameter space are expected to be probed directly by new exper¬ 
iments, and are already limited indirectly from astrophysical observations as shown in Figure 
Among the direct search experiments are the light-shining-through-wall experiment ALPS-II |48j . 
the helioscope lAXO [1^, and the observatories PIXIE [50] an PRISM |5l]. Indirect astrophysical 
limits excluding portions of the parameter space are obtained from massive stars [62], the 1987A 
supernova |63II65| , and quasar polarization |661 [67] . Since the present limit coming from the super¬ 
nova 1987A |65) is stronger than the limits coming from quasar polarization, we do not show the 
latter in Figure]^ 


D 

Va 

GeV] 

9a'y / Caj 

[GeV-1] 

7 

6.04 

X 

10® 

1.92 

X 

10- 

-12 

8 

2.44 

X 

lOlo 

4.76 

X 

10- 

-14 

9 

3.43 

X 

lOii 

3.39 

X 

10- 

-15 

10 

2.50 

X 

1012 

4.64 

X 

10- 

-16 

11 

1.18 

X 

1013 

9.87 

X 

10- 

-17 

12 

4.07 

X 

1013 

2.86 

X 

10- 

-17 

13 

1.12 

X 

1014 

1.03 

X 

10- 

-17 


TABLE I: 


Necessary values for Va to obtain rua = 7.1 keV in Eq. (18) through the operator in Eq. (16), 


with g = 1. 


Next we present a general analysis of the symmetries and find requirements for an acceptable 
model. After that, we show specific constructions implementing the ISS mechanism with scales 
originating from the vev of one or more scalar fields, which contain ALPs having values of {ga-y, w-a) 
in the regions that could explain certain astrophysical phenomena. 
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FIG. 2: Figure adapted from Ref. [20j. Values for the ALP coupling, Qa-y, and mass, rua, required to explain 
indications of the anomalous Universe transparency for very energetic y-ray |30j , the soft X-ray excess from 
the Coma cluster [351 [35], and the X-ray line at 3.5 keV [33113011551153], are inside the regions delimited by 
red lines. In green are the prospective regions to be reached by the light-shining-through-wall experiment 
ALPS-II [48], the helioscope lAXO [43], and the observatories PIXIE [50] an PRISM [51]. Also shown are 
astrophysical limits provided by: emission of ALP from massive stars representing an anomalous energy loss 
and shortening their helium-burning phase so that Cepheids could not be observed [53], non-observation 
of a y-ray burst emitted by the core of the supernova SN 1987A in coincidence with their neutrinos burst 
arrival on Earth [53H55] . The benchmark points for the one ALP model of Section [lIB| (A, square) and for 
the two-ALP models of Section [Tl C l| (B.l and B.2, stars) and Section [TlC 2| (C.l and C.2, triangles) are also 
shown. 


A. Symmetries of the model 

Two symmetries are essential in our construction: the continuous anomalous U(l)x symmetry 
and the stabilizing discrete gauge symmetry Z^v- Here we consider a single factor for simplicity but 
more factors can be equally considered. The continuous U(l)x symmetry should arise accidentally 
from the conservation of Z^r at the intermediate ALP scale. The discrete symmetry, in turn, is 
assumed to be a remnant of a continuous gauge symmetry at higher energy scales, possibly at the 
GUT scale [8] . We will be concerned with the theory at intermediate scales and we will not attempt 
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to trace the original continuous gauge symmetry as the possibilities are numerous. The anomalous 
nature of U(l)x gives rise to the required ALP-photon coupling whereas the discrete symmetry 
Z^r should be anomaly free in the discrete sense [U] as required by its local nature. 

Let us proceed to determine the generic aspects of the anomalous U(l)x and the discrete Z^r 


symmetries. We will establish the following: U(l)x acting on non-SM fields as in Eq. (12) is an 
extension of the usual lepton number L and Z^r can be chosen as a discrete subgroup of some 
combination of C and the baryon number B. 

We start with U(l)x and consider the Yukawa interactions of th SM: 


= ciLHdR + cilHur + LHIr , 


( 20 ) 


where we omit Yukawa couplings and family indices for simplicity. The SM helds are denoted as 
follows: qn are the left-handed quark doublets; UiR and diR are the right-handed quarks singlets; 
liR are the right-handed lepton singlets; with Li and H being respectively the left-handed lepton 
and Higgs doublets of Eq. ([^ . 


There are three family independent U(l) symmetries in Eq. (20), coming from 3 independent 
constraints on 6 phases associated to 6 types of fields. They can be identified as hypercharge Y, 
baryon number B and lepton number £. We adopt the usual assignment that the lepton fields 
Lj, liR carry one unit of lepton number: £ = 1. 

We now consider the addition of the right-handed neutrino fields Nr, Sr, necessary for the ISS 
mechanism, and also the complex singlet scalar a whose vev sets the neutrino mass scales M, n. 
These three complex fields contribute to the Lagrangian in Eq. 0: containing three terms, and no 
additional U(l) symmetry appear. These new fields do not carry neither hypercharge nor baryon 
number, the latter following from the absence of interactions with quarks. Thus they carry an 
extended lepton number. In particular, because of the first term in Eq. 0> NiR carries the same 
lepton number as Li. If we denote by a,d the lepton number of Sr and a, respectively, the last 
two terms in Eq. ([^ result in 


1 


d= {p- 2?) 


o = Iqd , 


( 21 ) 


where p ^ q/2 is required from independency of constraints. If we exchange a by a* in the last 
term in Eq. 0> it is sufficient to consider negative g —)• — q in all equations. We conventionally 
adopt positive p. 

At last, the new vector-like fermion fields Er,Er have the same electric charge as liR and its 
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hypercharge is defined. They interact through 


-££ D ki—^LiHEn + kE _^ ElEr , 

iWpi Mpj 


( 22 ) 


with small integers r, s; note that s cannot be zero to generate an anomalous symmetry but it 
should also obey |s| < 3 to induce sufficiently large masses for E, for ALP scales in the range in 
Eq. (§. Roughly speaking, the value of r determine the life time of the charged lepton E: the 
larger the value of r, the longer the life time of the E particle. If r is too large, the first term in 


Eq. (22) becomes negligible, and thus the E exotic lepton will be a stable charged particle which is 
cosmologically problematic, unless its mass is < TeV m- Another constraint comes from searches 
for long-lived charged particles in pp collisions m- We will discuss this in more detail below. 


The two interaction terms in Eq. (22) then determine the lepton numbers of Ei,Ee without 
affecting the number of symmetries. If we denote the lepton number oi El, Ee by b, c respectively, 
we obtain explicitly 


b — c = sd, 


c = 1 — rd. 


(23) 


A negative s,r in Eq. (23) may account for the simple exchange cr —)• a* in the respective terms. 


The final set of U(l) symmetries of the model consists of Y,B,C generated by charges listed in 
Table [11} The T-charges a,b,c,d of fields Sr, El, ER,a were determined by Eqs. (21) and (23). 
Since B, and obviously Y, are not anomalous with respect to electromagnetism, the anomalous 
symmetry U(l)x can be chosen to be generated by some multiple of the extended lepton number 


C. Specifically, since the anomaly is proportional to 6 — c, s cannot be zero in view of Eq. (23). 
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TABLE II: 


Continuous symmetries of the model. Charges a,b,c,d are determined by Eqs. (211 and (231. 


Concerning the discrete symmetry Z^r, the following anomaly cancellation conditions should 
hold from the effective point of view 13 do]: 


^2 = ^3 = ^grav = 0 mod N/2 . 


(24) 
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where A 2 , As, Agj-a_^ are the anomaly coefficients associated with [SU(2)i]^ x Ztv, [SU(3)c]^ x 
and [gravity]^ x Z^r, respectively. We ignore the anomaly associated to [U(l)y]^ x Z^v because it 
does not furnish useful low energy constraints [HI b] 

We write the action of Z^r as 

^ , (25) 

with discrete charges = Z{'il)k) = 0,1, • • • ,N — 1. Given that U(l)y is anomaly free by con¬ 
struction, and its imposition has no effect on undesirable operators, we can consider our discrete 
Zat to be a discrete subgroup of the rest of the symmetries in Table |n][l6l[l8]: 

Z = ciB + C 2 C , (26) 

where Cj should be rational numbers that makes all Z charges integers. To avoid redundancy, we 
can adopt ci = ni3 and C 2 = n 2 C 2 , where ni, n 2 = 1,..., — 1 and C 2 is the smallest integer that 

makes all £-charges integer and coprime. The factor 3 in ci appears because only 3B is made of 
integers. If N is not a prime, we also need to discard values for m that makes ciB or C 2 C to have 
a common factor that divides N. The latter case implies only a subgroup of Z^r is realized. 

We can now calculate the anomaly coefficients as 

Ai{Z) = csA{B) + C2Ai{C). (27) 


where 


AiB) = (1,0,0), 


(28) 


-^i(B) = (|) 0, —3a + 6 — c), 

with i = 2,3, grav, respectively. In special, the gravitational anomaly only depends on C and we 
can write 


^grav 


(Z) = C2(s - h)d. 


(29) 


where Eqs. (21) and (23) are used. We can see the well-known result that B — C is anomaly free for 
a = 0 and b = c, which corresponds to the SM with three right-handed neutrinos; see e.g. Ref. [72]. 
Therefore, any discrete subgroup of B — C will have ^42 and As automatically canceled. However, 


due to its discrete nature, Eq. (24), distinct combinations of B and £ can be also anomaly free as 
well. One can also check, there is no intrinsic discrete symmetry besides subgroups of combinations 
of Y, B, C; use, e.g., the Smith Normal Eorm method m- 

To summarize, we seek SM extensions defined by Eqs. 0 and ( |22[ ), with ALP decay constant 
Va, integers (p,q,r,s) and discrete symmetry Z^r obeying the following restrictions: 
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1. One ALP is present that couples to photons and explains one or more astrophysical phe¬ 
nomena indicated in Fig. 

2. The correct mass scales for the ISS mechanism are generated by 

3. The ISS mechanism is stabilized by Zjv; 

4. The mass of the heavy lepton E is larger than the electroweak scale: Me ^ r’w 

5. There is no discrete anomaly for Z^v- 

Extensions to more than one ALP should obey analogous conditions. 

The conditions for items 1, 2, 4 and 5 have already been discussed. To summarize conditions 
1 and 4, it is necessary to have 0 < |s| < 3 and the singlet a should be charged by U(l)x; 
and |s| = 3 is possible only if the ALP scale is high, > 10^^ GeV. The stability of the ISS 
mechanism, condition 3, requires the following: Nr, Sr should be charged under Z^r to avoid the 
direct Majorana terms NrN'^ and SrS^. Moreover, Z^v charges should prevent the appearance 
of operators a'^SnS^, a^SnN^, u'^NrN'^, u'^LHSr with dimension lower than the ones inducing 
the correct ISS scales in condition 2; and the same applies for operators that replaces a with 
a*. Specifically, any operator of the form a'^NRN^ or a'^'LHSR disrupts the zeros in the ISS 
texture in Eq. Q. Nevertheless, the mass matrix in Eq. Q is the leading contribution as long as 
W\ > \q\ “P) I^I > 2p— |g| and |n| + \n'\ > |g|; order of magnitude conditions can be extracted from 
subleading contributions in the seesaw formula, cf. [51 El]. To guarantee that these contribution 
are negligible, we require a more strict condition: \n\ > 4 and \n'\ > 3. The presence of these 
dangerous operators can be traced from their £-charges: 

LHSR^{q-p)d, N^RNR^{2p-q)d, (30) 


where we have conveniently written the charges in terms of the charge of a ~ d. Therefore, the 
combinations q — p and 2p — q control the coupling of these operators to powers of a and some 
combinations of p, q can be readily excluded in the case of one singlet. For example (p, q) = (2,3), 


is excluded because it allows both operators in Eq. (30) to couple to a*. Generically, it is more 
interesting to have negative q when p ^ 0- 


B. Model with one ALP 


We focus first on a model which according to our previous considerations could explain the 
Universe transparency and the soft X-ray excess from the Goma cluster. As we pointed out in Figure 
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there is an overlap in the parameter space and for certain values of {ga-y, ma) the same ALP could 
be responsible for both phenomena. In addition, to obtain correct order of magnitude parameters 


for the ISS, we take the singlet vev to be within the first interval in Eq. (10), corresponding to the 


case (p, q) = (2, 3). This choice leads to an ALP-photon coupling constant in Eq. (14) in the range 


9a'y 


2.1 X — 5.3 X 10“^^ GeV^ 


(31) 


a7 


For a coefficient of order one, the value of Qa-y would be outside the region required to explain 
the Universe transparency. But it would be still possible to explain the soft X-ray excess from the 
Coma cluster if Ca'y ~ 5 and the ALP mass is restricted to nia < 10“^^ eV [35l [36] . 

We choose the Lagrangian given by Eqs. 0 and ([^, with {p,q) = (2,-3), (r, s) = (1,2). For 
convenience, we write explicitly only the terms that depend on the singlet a: 


-C D 




1 


< 7*3 ^ 


Vijj^^SiRiNjuy + ^Cijj^SiniSjR) 


+ ki 


(y ^ 




—+ %—+ H.c 


(32) 


The choices above lead to d = 2/7 in Eq. (21) for the lepton number of a. As a consequence, only 
7£ has all charges integer and the stability of the ISS mechanism requires either a discrete symmetry 
Zii or Z 13 . Other choices allow operators of the forms LHSrct'^^N'^Nrct'^ with dimensions that 
are too low. However, because of gravitational anomaly, we choose Z 13 as it can be seen below. 
The simplest possibility for the anomaly free discrete Z 13 symmetry is 


Z = 6 H + 7£ , 


(33) 


where the coefficients in Eq. (26) are chosen as follows: C 2 = 7 is kept and ci = 6 is chosen from 
the cancellation of mixed [SU( 2 ) 2 ,]^ x Z 13 anomaly. The explicit anomaly coefficients for Z 13 are 


A 2 — 39/2, A 3 — 0, Agrav — 13, 


(34) 


which are all zero modulo N/2 = 13/2. Hence, since the gravitational anomaly only depends on £, 
we can see only Z 13 is anomaly free and we can discard Zn. The ISS mechanism is also stable as 
the lowest order Z 13 invariant operators that could disrupt the mechanism are g^LHSr. 


The explicit charges Zj for each field can be seen in Table HI Note that, for the SM model fields, 
the Zi 3 charges are equivalent to 6 (H — £) modulo 13. 


We can see that the imposition of the Z 13 symmetry in Eq. (32) successfully leads to an acci¬ 


dental U(l)x symmetry, corresponding to the extended lepton number £ in Table [TT] with charges 
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TABLE III: Z 13 charges Zi 


{6B + 7C)i in the notation of Eq. (251. 


conveniently rescaled in Table IV to give X^j = 1. Such a symmetry coincides with the usual lepton 
number for the SM fields, but it is anomalous for El, Eji fields. The accidental U(l)x symmetry 
is only approximately conserved because it is explicitly broken by higher dimensional Z13 invariant 
operators suppressed by the Planck scale. 
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TABLE IV: U(I)x charge assignments for the fields in Eq. (32). 


As the singlet field a acquires a large vev, Uo- ~ 3 x 10^'^GeV, the anomalous U(l)jv is spon¬ 
taneously broken, making its phase field a{x) in Eq. (§ a pseudo Nambu-Goldstone boson. The 
singlet vev will also generate the ISS mass parameters in Eq. ([^ in the correct order of magnitude. 
The characteristic shift symmetry for a{x) is broken by Z13 invariant operators of the form in Eq. 


(16), where is the operator of smallest dimension. The latter operator gives the dominant 


contribution in Eq. (18) for the ALP mass, with magnitude 

rua ~ 1.58 X 10“^®eV 2 ^ 


3 X IQio GeV 


5.5 


(35) 


With the charges in Table IV and from Eq. (15), we can readily calculate the anomaly coefficient 
Ca'y = 4 and the ALP-photon coupling, 

'3 X IQi^GeV' 


da-y 


1.5 X 10 


-13 




GeV 


-1 


(36) 


The benchmark point for this model, named A, is shown in Eigure Although the ALP in this 
model may explain the soft X-ray excess from the Goma cluster, it is out of the projected regions 
for searches of the ALP-II |l8] and lAXO |l9] experiments. But it is inside the region that will be 
probed by the planned observatories PIXIE/PRISM [5011^ . 

Since all U(l)x and Z13 charges are family blind, our model does not lead to specific predictions 
for the neutrino flavor structure (family dependent U(l) symmetries leading to axions has been 
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considered in, e.g., ref. m)- Only the order of magnitude for the absolute mass scale is obtained 
through Eq. 0 as 

m-i/~ [y^ x 1.4x lO^eV, (37) 

for Uw = 246 GeV and = 3 x 10^^ GeV. The mass matrices in Eq. Q for the ISS mechanism 
acquire quite natural values as 

M = r] X 187 GeV , /x = C x 1.6 keV . (38) 


Typically we will need the matrix entries of r] to be larger than one, e.g., ry ~ 5, and y <C 1, to 
evade lepton flavor violation [55] and direct detection constraints |56j as well as to maintain the 
validity of the seesaw formula in Eq. Q. On the other hand, C can be of order one or smaller. 
To obtain light neutrino masses in the sub-eV range, we need small Yukawa coupling, y < 0.1 or 


smaller. The largest eigenvalue of the combination of matrices inside brackets in Eq. (37) needs to 
be at most around 10“^. 


From Table IV we note that the charged lepton E is the only one that contributes to the 
coefficient Ca'y because E has electric charge different from zero and Xej^ ^ Xe^- In this model, 

2 _ I 

its mass comes from the term in Eq. (32). When a gains a vev, E obtains a 

2 

mass. Me = ^ kE x 187GeV. Therefore, we typically need kE to be larger than one to 

avoid the current lower limit of 574 GeV on the mass of new charged leptons m- To be more 
specific, the latter limit applies for charged long-lived heavy lepton with lifetime greater than a 
few nanoseconds, because these particles can travel distances comparable to the size of modern 
detectors and thus appear to be stable. However, in this model E can decay into and h^, with 
i = e, y, T and being the Higgs boson with mass of m/jO = 125 GeV. This decay is induced by 


the term kij^^LiHE^ in Eq. (32). Estimating the lifetime of E, te, we find that for ttle > 
te can be written as 


167r 


m'i 


te 


X 6.5822 X IQ-^'^s. 


3A2 

where we have neglected the masses of the SM leptons, i.e., me,m^,mT 


(39) 


0. We have also 


considered that ke = k^ = kr- The factor A in Eq. (39) is . Taking ke = k^ = k^ = 

find that for mE ^ 250 GeV, the charged lepton E has a lifetime smaller than 10“® s. Therefore, 
the lower limit of 574 GeV does not apply and order one values for kE are still allowed. 

Concerning other possibilities, a few remarks are in order. Considering the SM augmented by 
only one singlet scalar a and fermion fields Nji, Sr, El,Er through Eqs. 0 and ([22]), 
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• there is no other model that contains one ALP capable of explaining the transparency of the 
Universe or the soft X-ray excess in the Coma cluster; 

• it is also not possible to find models featuring an ALP with mass rua ~ 7.1 keV, which can 
explain the 3.55 keV X-ray line through the decay of the ALP into two photons |39t 140). 

Other choices for the powers {p,q), such as {p,q) = (3, ±5), and for the discrete symmetry Zjv do 
not comply with one or more of the restrictions explained in the end of Section p^I A[ (i) generation 
of correct mass scales for the ISS mechanism and (ii) stabilization of the ISS mechanism and 
(iii) cancellation of discrete anomalies. Many possibilities are excluded by (ii) because they allow 
low-dimensional operators a” to couple to N^Nji or LHSr. Further restriction comes from the 
gravitational anomaly cancellation, (iii), and only the Z13 symmetry model survives. 


C. Models with two ALPs 


Here we extend the previous setting and seek models featuring two ALPs that can explain the 
excess of X-ray photons in the 3.5 keV line, in addition to the transparency of the Universe for ultra 
energetic gamma rays and the soft X-ray excess from the Coma cluster. As we have previously 
discussed and shown in Figure at least two ALPs are necessary to explain these three phenomena. 
Hence, in addition to the singlet cr, we introduce another SM singlet a' which will host a second 
ALP a'. Now, two energy scales, Va = '/2{a) and Va' = \/2(cr'), will govern the physics of these 
ALPs. Since the ALP a' should also couple to photons, the singlet o' should be charged under 
another anomalous symmetry U(l)x', which follows accidentally from a second Zjv' symmetry. 

Let us choose a to be the ALP of 7.1 keV mass that explains the 3.55 keV X-ray line. The 
possible values for and the Z^r symmetry that are needed can be seen in Table [T] for \g\ = 1. The 
possible values for v„i are then restricted by Eq. (19), v^t Ri 10® GeV, and o' should be protected 


by a discrete symmetry Z^r/ with N' > 11. Note that the two ALP scales do not mix in our models, 
following the ones proposed in Ref m- Models where only o (or o') couples to Nr, Sr are excluded 
from the considerations of the previous section. We need that both and Uo-' generate the ISS 
mass scales. However, we were unable to find a plausible model that could satisfy all conditions 
listed in the end of Section p)! A[ Thus we present in the following, two models that satisfy almost 
all criteria. 

The general Lagrangian we will consider is composed of the usual Yukawa interactions in Eqs. 
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and 0, with terms depending on the scalar singlets modified to 


^Pl 


^Pl 


+ fcj 




jLiHEn + kE 


ctV'" 


( 40 ) 




■Ei^lEr + H.c. 


Many restrictions on the integers {p, q, r, s) discussed in Section IIA and in the beginning of Sec¬ 


tion |II] are now valid for the sum of unprimed and primed variables. For example, the restriction 
in Eq. ([^ should be now adapted to p -|- p' = 2,3 and |g| -|- \q'\ = 3,4, 5, where we conventionally 


take p,p' to be positive. Likewise, condition 4 in the end of Section IIA, for low ALP scales, is 
now |s| + I s'I < 2, which leads to s = ±1 and s' = ±1. We also see that the number of symmetries 
are consistent: there is one more field a' for the same number of constraints but we need one 
more anomalous symmetry. Given that the fields beyond the SM only couple to leptons, we can 
still consider X and X' proportional to two extended lepton numbers C and C!. Additionally, all 
formulas for (p, q, r, s) in Section [T^ still apply considering that a is only charged under C while 
a' is only charged under C. Thus the same formulas applies for the primed {p', q',r', s') as well, 
depending now on £'-charges (o', 6', c', d') of Sr, Er, Er, o' , respectively. 


1. Model I 


The first model gives up the cancellation of the gravitational anomaly for one of the discrete 
symmetries Zjv or Z^r/. The anomaly can be easily canceled by the addition of one or more fermions 
that are singlets of the SM but do not contribute to the phenomena discussed in this paper. The 
model also gives rise to a small scale for p. 

We consider the interaction terms for the singlet fields to be 


-.if D + CijJp^Ri^Rj 

. oo' — . oo' - 

+ ki , LiHEr kR-r-? — ERER-\-Yi.c. 


(41) 




'Mpi 


The Lagrangian has the form in Eq. (54) with {p,q) = (2,0) and {p',q') = (0,-3), whereas 
(r, s) = (r', s') = (1,1). We choose Va- ~ 2.44 x 10^*^ GeV to accommodate the correct scale for M 
and generate the ALP mass of rua = 7.1 keV, which should be protected by a Zg symmetry. We 
protect o' with a symmetry Zn so that the whole model has a symmetry Zg x Zn. 


Let us proceed to find the symmetry Zg x Zn. Eq. (21) determines {d,a) = (1/2,0) and 


{d',a') = (2/3, —1). Therefore C2 = 2 and c'2 = 3 makes all charges of C 2 C, and c'2£' integers and 


20 















we can calculate the gravitational anomaly from Eq. (29); 


^grav(c2>C) = 1, Agrav(c2^0 = 11- 


( 42 ) 


It is clear that the gravitational anomaly for does not cancel for any Z^v, except Z 2 . Note 
that we can not use C 2 = 2 x 4 to cancel the gravitational anomaly because 8/1 only generates 
Z 2 . Therefore, we assume such a gravitational anomaly is canceled by additional fermion fields 
and we adopt Zg generated by Z = 6;B + 2/2 (we could have adopted Z = — 6/2 as well). This 

choice cancels the anomaly of A 2 {Z). Analogously, we choose the Zn generator as Z' = —3i3 + 3/2'. 
We show the explicit charges in Table fV} The charges for the anomalous U(l)x and U(l)x' are 


presented in Table VI 
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TABLE V: 


Zg X Zii charges in the notation in Eq. (25). 
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TABLE VI: Charges for U(l)x and U(l)x'- 


This model yields Ca'y = Ca'^^ = 2, cf. Eq. (15), which leads to the desired ALP-photon 
couplings 


a 

2ttVc 

a 

2'KVrr' 


Ca'y ^ 

r 

s 9.52 X lO-^^GeV^ 

da '7 ^ 

s 2.32 X 10"^^ GeV^ 


(43) 


for Va ~ 2.44 X 10^° GeV and « 10® GeV. The ALP masses are given by 

~ \g\ X 7.1 keV, nia' ~ \g'\ x 3.41 x 10“^^ eV. (44) 

Benchmark points for ALPs a and a' are marked as B.l and B.2 in Pigurej^for \g\ = \g'\ = 1. 

The induced neutrino mass matrices have magnitude 

M = gx 124 GeV, g = Cx 0.061 eV, (45) 
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which leads to the light neutrino mass matrix 


m,. 


= [y V ^C(f? ) ^y] X 0.12 eV. 


(46) 


Analogously to the model of Section IIB we typically need rj to have entries with magnitude larger 
than one while the Yukawa coefficients need to be smaller than one, so that e = has small 

entries. The matrix ^ can have entries of order one or smaller but we can see the scale generated 


by a' is smaller than the one generated in Eq. (38). The ISS mechanism is stable as the new 
operators of lowest order are and LHSjia^a'^. 


2. Model II 


In the second model featuring two ALPs, we require the cancellation of all anomalies, including 
gravitational anomalies, but we relax the conditions for stability of the ISS mechanism. Because of 
the former, we can only find a symmetry Zg x Zio, so that the ALP a' is heavier than the previous 
model and can account for the y-ray transparency problem but not the soft X-ray from the Coma 
cluster. 

The model Lagrangian involving a, a' is 






+ h^UHEn + kE^EiER + H.c.. 


(47) 


'M|i 


'Mpi 


The Lagrangian has the form in Eq. (54) with (p,(?) = (2,-2) and {p\q') = (0,-1), whereas 
(r, s) = {r',s') = (1,1). As in the previous model, we choose Va ~ 2.44 x 10^*^ GeV and a is 
protected by Zg. The symmetry Zio then protects a'. 

After performing the calculations of Section II A| for this case, we choose Z = 3C — 3B and 
Z' = C! 9B as generators of Zg and Zio, respectively; they are given in Table VII[ One can 
check that these charges are anomaly free. The anomalous symmetries U(l)x and U(l)x' can be 
obtained from the same table by eliminating the baryon number contributions and rescaling the 
X' charge of a' to unity. The extended lepton numbers C and C! can be extracted in an analogous 
manner. They give Ca-,, = C'a/.y = 2. The ALP-photon couplings are the same as for model I, 


gar, ^ 9.52 X 10-1^ GeV-i , ga^r, 
for Va ~ 2.44 X 10^*^ GeV and v^' ~ 10® GeV. The ALP masses are given by 

\g\ X 7.1 keV, nia' ~ \g'\ x 1.81 x 10“^® eV. 


2.32 X 10"^® GeV^ 


rUn 


(48) 


(49) 
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We can see rua' is too large to explain the soft X-ray excess in the Coma cluster. Benchmark ALP 
photon couplings and masses can seen on Fig. [^marked as C.l and C.2. 
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TABLE VII: 


Zg X Zio charges in the notation in Eq. (25). 


Finally, the neutrino mass matrices have magnitude 

M = r] X 124 GeV, /r = C x 36.5 eV , 


(50) 


which leads to the light neutrino mass matrix 


= [y'^ v ^y] x72eV. 


(51) 


In this case, we have a more natural scale for fj, compared to the model of Section IIC 1 but the 
general considerations for rj are the same. Concerning the stability of the ISS mechanism, we can 
see the lowest order operators that disrupt the texture in Eq. (§ are and LHSRa*^a'. 

They lead respectively to mass parameters of magnitude 


y-N 


VgVa' 


36.5 eV, niDS 


VwVgVa' 
2^/2 M|i 


2.6 X 10"^^ eV. 


(52) 


These mass matrices contribute to entries (3,1) and (2,2) in Eq. (§ and contribute to subleading 
terms in the light neutrino mass matrix in Eq. © as isaEn 


— 

We can see the contribution of niDS is negligible. The mass parameter /ttv is of the order of y but 
it also contributes negligibly to the light neutrino mass matrix, even if one-loop corrections are 
taken into account m- 


III. CONCLUSIONS 

We have generically studied the construction of models where one ALP results from a scalar 
singlet carrying an anomalous extension of the lepton number L of the SM. The ALP successfully 
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accounts for some intriguing astrophysical phenomena such as the soft X-ray excess in the Coma 
cluster, and at the same time, the singlet vev furnishes the correct mass scales that implements 
the ISS mechanism for neutrino mass generation through gravity induced nonrenormalizable terms. 
Moreover, the approximate nature of anomalous £-number, the ALP mass and the ISS mechanism 
are protected from additional gravity induced terms through a gauge discrete symmetry. The 
additional beyond SM fields are minimal: three families of right-handed neutrino fields Ni^, SiR, 
one heavy singlet lepton E and one singlet scalar a are added. 

By requiring the stability of the ISS mechanism and the cancellation of the discrete gauge 
anomalies, only one model survives, and the discrete symmetry needs to be a Z13 subgroup of a 
combination of C and the baryon number B, the simplest being 6B -|- 7C. 

Simple extensions to models with two ALPs can be constructed by adding solely one more singlet 
scalar. In this case, two ALPs can solve more astrophysical phenomena with distinct features. We 
have been unable to find a model capable of explaining the three astrophysical phenomena and, 
at the same time, satisfying all the conditions in Section |il A[ Therefore, two models are presented 
by relaxing some of the conditions. The first model can explain all the astrophysical phenomena, 
but the gravitational anomaly for one Zjv factor can not be canceled within the field content, and 
additional fermionic fields are required. The second model does not present discrete anomaly but 
it can not explain the soft X-ray excess in the Coma cluster. In both cases, the 7.1 keV ALP may 
be an appreciable component of dark matter as well |20l [76] . 

In summary, we have proposed a very restrictive and economical setting to extend the SM and 
explain notable astrophysical phenomena together with natural neutrino mass generation through 
the ISS mechanism. The restrictiveness of the setting allows only one model with one ALP and 
models with two ALPs are also largely restricted. More possibilities emerge if we allow the presence 
of more than one heavy charged lepton and, in particular, three copies of them can lead to easier 
cancellation of the gravitational anomaly. 

Erratum 

The two models containing two ALPs in Sec. 2.C.1 lead to too light non-SM charged leptons E of only 
a few GeV. (The single ALP model of Sec. 2.C.2 does not suffer from this problem.) 

This problem can be amended without significant modifications by considering two heavy vector-like 
fermions E and E' instead of one. Instead of the last two terms of Eq. (40), (41) and (47), we should 
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consider the four terms 


-^Dk—L^HEn- 

M2j 




+ kEcrELEn + k'^a'E'j E'„ + H.c., 


(54) 


where we already fixed the r, r' powers without affecting any formula. One should use s = s' = 1 if needed. 
Now the heavy leptons E, E' have intermediate scale masses of order 10®'^^°GeV. The charge of the discrete 
symmetries are only modified for El.r and ^ in a predictable manner; see table below. Appropriate 
normalization of these charges leads to the PQ charges X,X' and the extended lepton numbers L,L'. Note 
that U{l)x \U{^)x'] is vectorial for E' (E) and that the first two terms of the Lagrangian above implies 
Er and E'l^ have equal charges {X, X', Z, Z', L, L'). The rest of the formulas and the phenomenological 
consequences remain unchanged. 



El 

Er 

E'l 

E'r 

Zs 

2 

1 

1 

1 

Zii 

1 

1 

3 

1 



El 

Er 

E'l 

E'r 

Zs 

3 

2 

2 

2 

Zio 

-1 

-1 

1 

-1 


TABLE VIII: Corrected charges for model I (left) and model II (right). 

Additionally, the lowest order operator coupling L-Sr for model II is LHSRa'^a' instead of LHSRa'^a'*. 
This implies mos is negligible in Eq. (52). 
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